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Abstract. Using non-equilibrium renormalized perturbation theory, we calculate the 
retarded and lesser self energies, the spectral density p{uj) near the Fermi energy, and 
the conductance G through a quantum dot as a function of a small bias voltage V, 
in the general case of electron-hole asymmetry and intermediate valence. The linear 
terms in uj and V are given exactly in terms of thermodynamic quantities. When the 
energy necessary to add the first electron (Ed) and the second one {Ed + U) in the 
quantum dot arc not symmetrically placed around the Fermi level, G has a linear term 
in V if in addition either the voltage drop or the coupling to the leads is not symmetric. 
The effects of temperature are discussed. The results simplify for a symmetric voltage 
drop, a situation usual in experiment. 
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Figure 1. (Color online) Scheme of the system. A QD (chxle) is connected to the left 
(L) and right {R) conducting leads, with a difference eV between the corresponding 
chemical potentials. Each lead hybridizes with the dot level providing a contribution 
/3yA to the resonant width A. 



1. Introduction 

Recently, there has been great interest in the conductance through one quantum dot 
(QD) for low applied bias voltage and temperature, searching for universal scaling 
properties [HEllSllllElEllT]. In experiments [H El [TJ E] a semiconductor QD or a 
single molecule is attached to two conducting leads (left L and right R), a bias voltage 
V is applied between the leads (see Fig. [1]), and the conductance G = dl/dV, where / 
is the current, is measured. For V = 0, the temperature dependence of the conductance 
has been found to be very well described by the same universal function G{T/Tk), where 
Tk is the Kondo temperature, even for systems with very different [IllH]- This scaling 
law had been obtained from numerical renormalization group (NRG) calculations of the 
impurity Anderson model [9] in the Kondo regime {—Ed ^ A and Ed + U ^ A, where 
the Fermi energy is set as 0, -E^ is the on-site energy, U the Coulomb repulsion and A 
the resonant level width). 

In the non equilibrium situation 7^ 0, the problem is much tougher theoretically. 
In this situation, exact Bethe ansatz results are available only for a simpler problem (the 
interacting resonance level model) [in] while NRG methods have serious limitations [H] . 
In some works, the spectral density Paiuj) calculated with NRG at equilibrium is used 
to calculate non-equilibrium properties, assuming that Pa{oj) is not affected very much 
by an apphed bias voltage V [121 ESI [IH [15] • However, this approach misses the effects 
of broadening of Pa{oj) caused by V , and as a consequence, the results are quantitatively 
and in some cases even qualitatively different than the correct ones [T6] . 

In the Kondo limit, for either eV ^ kT^-, or /i^-B ^ kT^ in presence of an applied 
magnetic field S, have been determined using perturbation theory and poor man's 
scaling [17]. Using a Fermi liquid approach, based on perturbation theory in U (PTU 
[T8l fT9l [20I [2Tj). and Ward identities, Oguri has determined exactly the scaling up to 
second order in kT and eV for the symmetric Anderson model [221 123] . Further work 
considered the effect of higher order contributions using different approximations, like 
PTU [3], non-crossing approximation [5], or decoupling of equations of motion [6]. 
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The effect of asymmetric coupling to the leads A/j [defined by Eq. ([2])], and 

asymmetric drop in the bias voltage [ai 7^ ctR, see Eq. ([2])] has been calculated up to 
second order in T and V using Fermi liquid approaches, for the symmetric Anderson 
model {—Efi = + U) [3l H]. In other works electron-hole asymmetry was present for 
V = Q but symmetric voltage drop and coupling to the leads was assumed [5l E]. In all 
these works, the linear terms in T and V are absent, and the conductance is maximum 
for T = y = 0. In this work, we calculate exactly the linear term in V in terms of 
thermodynamic quantities at equilibrium. This term is non vanishing if —Ed ^ Ed + U 
and in addition either 7^ A/? or 7^ oir- Previous work for the infinite U Anderson 
model, a expansion pi] predicts that the conductance has a linear term in V of 
the form 2(Ar - Ai)ey/(AfcTi^) for = ur. 

A comparison of recent theoretical studies of the scaling properties [HI HI [5] with 
experiment indicates that some degree of intermediate valence (a deviation from the 
Kondo regime) is necessary to explain the experimental data of Grobis et al. [1], 
in which the effect of V to decrease the conductance is about 2/3 smaller than that 
expected in the Kondo regime. A good agreement with this experiment is obtained 
using the Anderson model [3], [5], while the experimental results of Scott et al, [2] seem 
to be out of the reach of the model (see SectionH]). For semiconductor QD's, it is believed 
that the Coulomb repulsion U is not too large in comparison with the resonant level 
width A = A/j + Ai, [25]. In a recent experiment in which the shot noise is measured, 
these parameters were estimated as U = 0.56 and A = 0.17 meV respectively [26] . 

The above facts suggest that the expansion parameter in PTU U/ttA is of order 
one, and then PTU might be a useful approach. In fact, PTU up to second order in 
U (using an interpolative scheme to recover the atomic limit [271 EH]) has been shown 
to describe well the equilibrium conductance for fZ/vrA < 2 [29]. The results agree 
with those obtained using the finite temperature density matrix renormalization group 
method [30]. Interpolative PTU has also been used to calculate persistent currents in 
rings with an embedded QD or side coupled to a QD [31], giving results in agreement 
with exact calculations in small systems [31] , and renormalization group results [32] for 
U = 6.25A. Other applications at equilibrium include impurities in quantum corrals 
[33l[34], magnetoconductance [35], and conductance through Hubbard chains [36] . 

The extension of PTU up to second order to the non equilibrium (finite V) case 
has been first considered by Hershfield et al, [37]. They found that for finite V, the 
current is conserved only in the electron-hole symmetric Anderson model. Different self- 
consistent approaches were proposed to overcome this shortcoming, by a suitable election 
of the unperturbed Hamiltonian [27l [38] . While these approaches work well in absence 
of a magnetic field B, numerical difficulties were found for small non- vanishing B and V 
[38] . Applications of PTU up to forth order (which seem necessary to obtain a splitting 
of the spectral density for large enough V) [391 SO] , or calculations of the noise spectrum 
[m II2] were restricted to the symmetric model. We show in this work that the current 
is conserved up to terms of order without the need of adjusting the unperturbed 
Hamiltonian. We use renormalized PTU (RPTU) to calculate the conductance for low 



Nonequilibrium conductance of a nanodevice for small bias voltage 



4 



V. The basic idea of RPTU as developed by Hewson [13] is to reorganize the PTU in 
terms of fully dressed quasiparticles in a Fermi liquid picture. The parameters of the 
original model are renormalized and their values can be calculated exactly from Bethe 
ansatz results [H SSI SSI SH HH] , or accurately using NRG [l9l EQl pi |5| |53] . One of 
the main advantages is that the renormalized expansion parameter [//(vrA) ^ 1, being 
1 in the extreme Kondo regime {U — > oo). Here we assume B = 0. Exact results for 
small magnetic field were obtained using RPTU in the symmetric case [52], while for 
large B, calculations with interpolative PTU were presented [38] . 

We calculate the self energies and the spectral density near the Fermi energy for 
small values of the frequency u and bias voltage V, and the conductance for small V. 
The linear terms in u and V are given exactly in terms of the occupation at the dot, 
magnetic susceptibility and specific heat at equilibrium. The lesser self energy is given 
exactly up to quadratic terms. We also show that the current is conserved up to terms 
of order in PTU. 

The paper is organized as follows. In Section II we describe the system and the 
impurity Anderson model used to represent it. We also review briefly the formalism 
of the perturbation theory and the idea of RPTU. Section III contains the results of 
the calculations, and several limiting cases of interest in which the general expressions 
become simpler. Section SI contains a summary and a discussion. 

2. Model and formalism 

2.1. Model 

A scheme of the model is displayed in Fig. [H The QD interacting with two conducting 
leads is described by the spin 1/2 Anderson model. In general, to use PTU, it is 
convenient to split the Hamiltonian into a noninteracting part Hq and a perturbation 
H' as 

H =Ho + H', 

kva (T 
kvcF 

H' =^ {Ed - elfj) + U Ud^Udi, (1) 

where Uda = d^^d^, and u = L,R refers to the left and right leads, with chemical 
potentials 

fiL = aieV, fiB. = -aneV. (2) 

respectively, with + = 1. Similarly, the couplings to the leads assumed 
independent of frequency are expressed in terms of the total resonant level width 
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A = Ai + as 

A. = 7rY,\VkuMu - Eku) = /3.A. (3) 

k 

In general e"^^ is determined selfconsistently, except for the electron-hole symmetric 
case with 5 = 0, for which g^^^ = p7l[38]. 

2.2. Green's functions and self energies 

The one-particle properties of the system, including the current, are determined by three 
types of independent one-particle Green's functions, the retarded G^{uj), the advanced 
C^iu) which are the complex conjugate of and the lesser ones G'^{uj). The 

retarded Green's function of the electrons at the dot for spin a, can be written as 

""'^(^'^ = c.-.;„ + !a-S;(c.) ' 

In PTU up to second order in U, the retarded self energy is approximated as [38] 

Ki^) = - + UM + El'iu), (5) 

where contains the contribution of order U'^. 

The lesser Green's function can be written in the form [38] 

Gliu^) = IGLMP - S.^m) ' (6) 

where g^^{uj), g^^{oj) are the retarded and lesser noninteracting Green's functions, and 
S^(ci;) is the lesser self energy. 

For large values of ?7/A, ordinary PTU in U is not reliable and Hewson proposed 
to reorganize the perturbation series in terms of renormalized parameters (which we 
denote with a tilde) such that U /(vrA) ^ 1. Within RPTU, the low frequency part of 
G'^^^ioj) is approximated as [13] 

G'-M ^ % ~ , (7) 

where 



^ = [1 - dK/duj]-\ el; J = z[e:^f + S;(0)], 

A = ^A, S;M = ^SrM, (8) 

and the remainder retarded self-energy is defined as 

Sr (^) = - KiO) - iodK/dio. (9) 

In Eqs. (IZD and (El), Sj;(0) and dEl/doj are evaluated at u = T = V = 0. 

In Fig. [2] we compare the spectral density Po-(w) = —lmG^^^{u)/n evaluated within 
PTU up to second order with the corresponding noninteracting renormalized result 
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Figure 2. (Color online) Full line: spectral density of the Anderson model for U = 7 A 
and Ed = —U/2. Dashed line: corresponding result for non-interacting quasiparticles 
(see text). 

{zpq{u) = —lmG^^^{uj)/n with So-(w) = 0) obtained from Eqs. ([7]) and ([8]) for the 
symmetric Anderson model {£eff = 0)- We obtain z = 0.2744 for the parameters of 
the figure. This comparison has been made by Rejec and Ramsak for other parameters 
[55]. One can see that aheady non-interacting quasiparticles reproduce rather well 
the low-energy behavior. A similar comparison, for a case with non-trivial frequency 
dependent A(a;) has also been made [56]. showing also a very good agreement for low 
An accurate calculation of including high orders in the perturbation expansion, 

should lead in turn to a precise description of the spectral density at high energies, 
since nothing is lost in the reorganization of the perturbative expansion. However, 
the calculation of higher order terms becomes involved due to the need of considering 
counterterms and corrections to the parameters [43] . 

The spectral density evaluated at the Fermi energy is related with the occupancy 
by the Friedel sum rule [51] 

sin^(7r(nrf^)) 

Pa(0) = T , (10) 

which allows one to relate the effective dot level with its occupancy, using Eq. ([7]) 

£:^^ = Acot(7r(nrf,)). (11) 

^rem^^^ can be calculated perturbatively using a renormalized interaction U which 
is given by the vertex function F at the Fermi energy [l3l HQ] . In absence of magnetic 
field, the subscript a can be dropped. The linear term in the specific heat and the 
impurity contribution to the magnetic susceptibility (determined using Ward identities) 
are given by [13] 



7c = 27r2A;po(0)/3, 

X =((7/iB)2(l + ^po(0))/2, 



(12) 
(13) 
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where 

PoM = ^ (14) 

is the density of free renormahzed quasiparticles. Eqs. (fTU]) . ( IT^ . (fT5i) and an accurate 
knowledge of thermodynamic quantities from Bethe ansatz or NRG, permits a precise 
determination oi z = p(0)/po(0), and the renormahzed interaction through the Wilson 
ratio R = l + Upo{0). 

To order f/^, the renormahzed self energies can be written as [38] 

= U^ f de,de,de,- ^o(^i)Po(^2)po(.3) 



uj + e-i - ei - €2 + If] 

x[(l-/(6l))(l-/>2))/>3) 

+ />i)/>2)(l-/>3))], (15) 



= - 2z7rf7^ J deide2Po{ei)po{e2)po{ei + 62 - uj) 

x/(ei)/>2)(l-/(ei + e2-u;)), (16) 

where 77 is a positive infinitesimal and f{u) = Ylv I^A'^)f{'^ " l^v) is an average of the 
Fermi functions at the two leads [f{co) = l/{e^^'^'^ + 1)]. 

2.3. The current 

Using the Keldysh formalism EH] , it has been shown that the current flowing between 
the left lead and the dot is 

HcAl 
^' = ^ 

Similarly, the current with spin a flowing between the dot and the right lead is 

Aie^R 



j du {2if{u - /iL)ImGS(a;) + G<{u)] . (17) 



h 



j du [2if{uj - PR)lmG'-^{uj) + G<{uj)] . (18) 



Of course, since the current is conserved one should have II = Ir = I- Then, from a 
weighted average of both expressions, G"^ can be eliminated, giving j58] 

, 87r/3i/3HeA 



h 



dujp{uj)[f{uj - Pl) - f{uJ - IJ'r)]- (19) 
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3. Self energies, spectral density and conductance for small bias voltage 
3.1. Self energies 

It is easy to see that to evaluate the differential conductance G = dl/dV at zero 
temperature to order V^, it is enough to calculate the spectral density and self energies 
up to total order n in a; and V (all terms of the form uj^V^~\ with < m < /, / < n), 
since the effective interval of integration in Eq. (fT9l) is eV , because the Fermi functions 
f{<jj) become step functions 6{—uj) at T = 0. For the same reason, the intervals of 
the two integrations in Eq. f|T6|l are of order eV, and therefore one obtains S^(a;) 
up to terms of second order, taking po{e) at the Fermi energy e = 0. Performing the 
integration we find 



S<(u;)= -t7r[M0)?U' 

X {OiaLcV - uj)I3l{I31 + 2l3l){aLeV - uf 
+ ei-aneV - u)Pr{2PI + Pl){aReV + ujf 
+ ^[(1 + aL)eV - uj]I31I3r[{1 + aL)eV - uf 
+ ^[-(1 + aR)eV - uj](31(3r[{1 + aR)eV + uf. 

(20) 

Using Im(cj + iri)~^ = — Z7r5(a;) and Eq. (JTSl) it is easy to realize that the above 
arguments also work for lmZl(a;) and performing the two remaining integrations [after 
the use of the delta function 6{uj + 63 — ei — €2)] one obtains up to terms of second order 
in u and V 



ImE^(a;) = - -[po{0)fuV - 27weV + S{eV)], 
7 = - aRf3R, 

S = 7' + 3/3L/3ij, (21) 

as previously found by Oguri [23] . 

Note that higher order terms in the perturbation series in U would lead to terms 
of higher order in V. Therefore, the above results are exact. They can be used to 
test approximations. For example, from Eqs. fl20l) and fl2Tl) . one realizes that the 
approximation for the lesser self energy S^(u;) ~ 2i/(a;)lmS'"(u;) used in decoupling 
approximations [591 El [60], is only valid for V = (or U = 0). 

From Eqs. ([6]), ([7]), ([H]), (HZD and ([IH]) one obtains 

1^-1^ = -^ J duj\G2io)\'[2f{oo)lm^'{io)+tt<{uj)]. (22) 

The integral introduces another factor eV. Thus, in order to obtain II — Ir to order 
V^^, one can replace [(^^(a;)^ by its value at w = and move it outside the integral. 
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Using Eqs. ( 12U]) and the integral turns out to vanish. Then 1^ — In = OiV"^) and 
the current is conserved up to order V^. 

In contrast to and ItoE^^u), an accurate calculation of the real part of the 

retarded self energy is more difficult. As it is apparent from Eq. (fT5|) or alternative 
expressions [371 138], the calculation of ReS'^'(a;) involves an integration at high energies, 
for which the approximation given by noninteracting quasiparticles [Eq. ([7]) neglecting 
is not good enough (see Fig. [2]). Fortunately S''(0) = and dTi'' {u) / du = 
by construction [see Eqs. ([H]) and ([9])], while dT,^/deV can be calculated from Ward 
identities: using the results of Oguri [23], one sees that dT,^ /deV = 'jdT,'^/dfi, where /i 
is a shift of both chemical potentials. In addition, using Eqs. (3.15) to (3.18) of Ref. 
[19], and dY7 /du = 0, one obtains dY7 /dp = Pq{0)U. Combining both results one has 



previous section. Note that for {uda) = 1/2 and U — )■ +oo, corresponding to the Kondo 
limit, X = 1 [l3]. This means that in this limit, the position of the Kondo resonance at 
equilibrium remains at the Fermi energy [Eqs. (I7j), ( fTTj) . S^(0) = and dY7' /dp = x] if 
the chemical potential is shifted, as expected from the common wisdom on the Kondo 
resonance. 

3.2. Spectral density and conductance 

Using Eqs. ([7j), ffTOl) . f|TT]) . fl2T|) and (123|) one obtains up to quadratic terms in u and 
eV, calling for brevity c = cos{n{niia)) and s = sin^n {uda)) 



C =s2[4c'-l + x2(l-2s2)/2] 

- scAd^Rei:''/du\ 

D =^xs^[2-8xc^ -x{l-2s^)] 

- 2scAa^ReS79w9ey, 

E =x^s^[-f\Ac' -l) + 6{l-2s^)/2] 

- scAd^ReE^'/dieVY. (24) 



This equation might be regarded as extension of the Friedel sum rule Eq. ffTOj) to finite 
small frequency and bias voltage. The annoying terms in ReS^ enter only the second 
order terms and vanish for the symmetric Anderson model. Inserting the above result 
in Eq. (fT9|) and deriving with respect to V the conductance up to terms of order is 
obtained 



9S79er = 7X. (23) 
where we call x = Pq{0)U = R — 1, where R is the Wilson ratio mentioned in the 




A 



2sc, B = — 27XSC, 
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G{V) = G(0) 



' A 



M 



V A / 



G(0) 

L 

M 



2sc{aL - oiR- 27a:), 
C{al + a\) + ^D{aL 



(25) 

(26) 
(27) 
(28) 



This is the main result of this work. Except for the second derivatives of ReS*", the 
remaining quantities are given exactly in terms of the characteristic energy scale A, the 
total occupation n = 2{naa) at the dot, the Wilson ratio R = x + 1, the ratio of the 
couplings to the leads f3i, and the distribution of the potential decays a^. The energy 
A can be obtained from thermodynamic quantities. For example from Eqs. ( ITTl) . ( |T2i) 
and (fM]): 

^sin2(7m/2). (29) 

37c 



A 



Up to second order in the interaction, d'^KeH'^ / du"^ has been calculated by Horvatic 
and Zlatic [201 El]. However, for the reasons discussed above, we do not expect this 
(rather complicated) result to be accurate enough for strong interaction. 



3.3. Particular cases 

The total occupation of the dot n can be controlled by the gate voltage Vg. Assuming 
that the coupling to leads is not affected by Vg, the occupation n = 1 maximizes G{0). 
However, experimentally, the scaling properties of G{V,T) have been studied not only 
for Vg that maximizes G(0), but for other values as well [1]. If one starts from the 
symmetric Anderson model, then by symmetry n = 1 for V = 0, c = cos(7m/2) vanishes 
and with it all terms involving second derivatives of ReS^. In this case, Eq. (|25|) 
simplifies to 



G{V)/G{0) = 1 - ^[(2 + x'){al + 4) 
2A^ 

- 37x(2 + x){aL - an) 

+ 9x\^' + f3M], (30) 

which is equivalent to the result found previously by Sela and Malecki using a different 
approach [1]. 

To model the voltage drops, in some works [3l [25] it has been assumed that the 
shift in the average chemical potential with applied voltage A/i = 'jeV = [see Eq. 
(I2TI)]. However, even in cases with very asymmetric couplings, the observed structure 
of the diamonds in G(y) indicates that ai ~ aji [611 [62]. Taking = aji, then 
7 = {/3l — Pr)/2, and the coefficients of Eq. (l25l) simplify to 
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- scA[d^Ref:/duj'' + 3(9'ReS79(e\/)']. (31) 

Note that in a 1/A^ expansion of the infinite U Anderson model, the value L = 
-2{Pl - Pr) has been found [21]. 

3.4- The maximum of G{V) 

In the general case for which the linear term in V, M ^ 0, the maximum of G{V) 
does not lie at V = 0, but at Kmx = -L/{2M), with Gmax/G(0) = 1 - LV(4M). We 
discuss the case ol = o.r and 2/3<n<4/3for which the curvature of the unperturbed 
quasiparticle density of states is positive [see Eqs. f lTT]) and (JH])]. In this case, because 
L ~ X, and M < 0, the largest possible Gmax is obtained for high interaction x = po{0)U. 
To estimate an upper bound for Gmax, we neglect the second derivatives of ReS*" in Eqs. 
(!3T|) and take the maximum possible value of U /vrA, namely one, which implies x = s"^ 
[see Eqs. ([II]) and ([H])]. This gives 

G^ax < 4Ps — [1 + i6,4_7,2_3 + 24P.vJ- ^^^^ 

where the product P = This function always increases with increasing P. Since 

the maximum value of P is 1/4, one has more simply 

G^ax<^sin2(7m/2), (33) 
for ai = OiR and — 1| < 1/3. 

3. 5. Effect of temperature 

For n = 1, the effects of temperature T and voltage V on the spectral density and the 
conductance are additive up to total second order in u, T, and V. The corrections 
to the self energy and conductance for = in this case, were discussed previously 
P H [221 [23]. In particular 

p(a;,T)_^ {2 + x')u:' + {TTxkTl ^^^^ 



p(0,0) 2A2 

l + 2x^). (35) 



G{T) ^ (TikT) , . 



G{0) 3A2 

For n 7^ 1, the corrections to the self energy in second order in U and up to second 
order in u and T were calculated by Horvatic and Zlatic [20l [21]. There is however in 
general (for V ^ 0) a term linear in T which comes from the term of order U in the 
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self energy [see Eq. (^], or 9ReS''/9T = Ud{nj;^) /dT. An accurate determination of 
the last derivative is not possible with the knowledge of the Green functions only for 
frequencies near the Fermi energy. For = 0, d{nda)/dT = because the system is a 
Fermi liquid and the lowest order correction to the occupancy goes as [H]. 

Including only terms linear in u, V and T, and using the results of the previous 
section one obtains 

p(uj) ^ sinfvm) ^ zU d-n? ^^^^ 

G(T) sinfvrn) , 

= 1 + ^^[eV{aL -an- 2^x) 

where the derivative is evaluated at V = T = 0. 



4. Summary and discussion 

Using Fermi liquid properties and non equilibrium renormalized perturbation theory, 
we have calculated the self energies and spectral density near the Fermi energy, and 
the conductance for small bias voltage V compared to the characteristic energy scale 
A. We have determined exactly the linear term in V in the conductance G, in terms 
of the A, the occupation of the dot n and the Wilson ratio R = x + 1. This term is 
different from zero if 7^ 1 and if in addition either the voltage drop is asymmetric 
(oi 7^ an) or the coupling to the leads is asymmetric {/3l 7^ /3r). It is important to 
interpret experiments, because due to the measurement method, there is a small voltage 
offset [61]. For ai = aji, the sign of this term and its dependence on Pl — Pr agrees with 
previous results based on a 1/A^ expansion of the infinite U Anderson model [21]. In 
addition, the observed GiV) is asymmetric with a term [61], which is unfortunately 
beyond the validity of our approach, because the current is conserved only to order V^. 

For an accurate calculation of the terms quadratic in V when n 7^ 1, it is necessary 
to know second derivatives of the real part of the self energy with respect to V and 
frequency u. These might be calculated combining numerical-renormalization- group 
(NRG) calculations [53] with renormalized perturbation theory [231 US]- For n ^ 1 
and finite bias voltage V , there is a correction of the real part of the self energy with 
temperature T, proportional to dn/dT, which is also beyond the reach of low-energy 
expansions. For ^ = 0, dn/dT = [4^, but this is not necessarily the case out of 
equilibrium. 

Values of the Wilson ratio R for the asymmetric Anderson model can be found for 
example in Table VIII of Ref. [51]. In the extreme Kondo limit, i? — ?■ 2 (x — 1) and 
n — )■ 1 simultaneously. Then, the effects of asymmetry and the coefficients of the second 
derivatives of the real part of the self energy tend to vanish. 

The shape of the diamonds in experiments suggest a symmetric voltage drop 
{a^ = aji) [62], except in arrangements like that of a scanning tunneling microscope, for 
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which there is a very asymmetric couphng of the nanoscopic systems to the conductors 
{Pl/Pr ^ 1 or ^ 1) . In the latter case, the conductance is very small, due 

to the factor PlPr in the expression for the conductance [see Eq. (l25l) and ( 126|) ] . For 
o^L = oiR and |ra — 1| < 1/3, we obtain the maximum of the conductance remains below 
the quantum of conductance times sin^(7r?7,/2), neglecting the effect of second derivatives 
of the real part of the self energy. 

These terms disappear for n = 1. This situation can be searched experimentally 
adjusting the gate voltage in order to obtain the maximum equilibrium conductance 
(V^ — )■ 0) . In this limit, our results coincide with those obtained previously, using a 
different approach [1]. If in addition, the voltage drop is symmetric, the conductance 
up to second order in T and V has a very simple expression |63] 



G = 4/3,/3^^[l-i^(l + 2x^) 

2 



-^^(l + 5x^)]. (38) 

The ratio of the coefficient of {eVf to that of {kTf lies between 2,/{2ti'^) = 0.152 
in the strong coupling limit (x = 1) to half of this value in the non-interacting case 
(x = 0). Reported values are 0.10 ±0.015 (Ref. [1]), 0.15 [7] and 0.051 ±0.01 (Ref. [2]). 
The latter seems inconsistent with the predictions of the impurity Anderson model. 

Note added: After acceptance of this work we became aware of Ref. [6l], which 
addresses a similar problem for the particular case 7 = (a^ = using a perturbative 
approach in £eff- According to the results presented, this approach up to terms of 
total second order, leads to the same equations as in the Ng approximation [591 EQ] 
[T,^{uj) = 2if{uj)lmS'^{u)], which trivially leads to the conservation of the current [see 
Eq. (I22l)]. but is unfortunately incorrect, as discussed in Section \3A\ 
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